HARD QCD: SOME RESULTS 
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I recall some points on the present status of QCD results in the short distance 
regions and illustrate the case of event shape in QCD radiation 



1 Hard QCD studies 

QCD studies have produced successful results, specially for observables in- 
volving short distance regions, for over a quarter of century. High accuracy 
calculations of hard distributions are possible and often one is able to make 
quantitative predictions which are strongly constrained. Calculations are of- 
ten difficult but there are constant improvements in techniques and this is 
always associated to a better understanding of QCD features. A consequence 
of these successful results is that QCD is often considered an established and 
understood theory. Actually, there are a variety of unsolved theoretical and 
experimental questions in QCD. They range from the two extreme regions. 
On one side we still need to understand how quark and gluon elementary fields 
change into hadrons at large distances. Since we think we almost dominate 
the short distance regions by perturbative (PT) methods, there are attempts 
to extend these methods and even the language to the study of non per- 
turbative (NP) corrections. This is an approach to large distances from the 
"short distance window" . Going in the other extreme region, the standard 
understanding is that in the ultraviolet region the theory becomes free and 
no structures are present. But, again, this is based on PT studies. Actually, 
it has been argued that short distance structures could be present in QCD u. 

Hard QCD is systematically and analytically studied by resummation of 
the PT expansion. This procedure starts from free quarks and gluons and 
reconstructs their interaction by using the point-like vertices in the QCD 
Lagrangian. It is reassuring to discover afterwards that quarks and gluons 
cannot be free at large distances, although this region it is not accessible to 
PT methods. In order to resumm PT expansions for physical observables 
one needs in general to know multi-parton distributions. Of course we need 
approximations. I recall that the present PT status of the art for multi- 
parton distributions concerns: collinear singularities, soft singularities, fixed 
order exact matrix elements. 

1. Collinear singularities. In general, all leading and next to leading 



contributions to the multi-parton .distributions are obtained by generalizing 
of the DGLAP evolution-equation □ (for a generalization to "exclusive" multi- 
parton distributions see a). They are given in terms of the anomalous dimen- 
sions which are known to leading and next-to-leading order. Recently there 
are attempts to compute NNNL contributions a. Singular contributions to 
multi-parton distributions require a collinear cutoff. Universality of collincar 
singularities allows a unified description of hard distributions for all processes 
such as lepton-lepton, lepton-hadron and hadron-hadron collisions. 

2. Soft singularities. Also in this case all leading (and some next-to- 
leading) contributions to the multi-parton distributions are known. They are 
given in terms of factorized soft radiation from hard emitters u. In the case of 
two or three emitters it is possible to factorize the soft radiation at two-loop 
U. For more emitters the two soft parton factorization requires the planar 
approximation (leading order in l/N c ). I recall the most important results 
and features of QCD we learned in the 80th and 90th: 

• coherence of soft radiation; 

• colour connection and angular ordering; 

• gluon multiplication and hadron multiplicity; 

• hump-backed plateau; 

• jet-finding algorithms in e + e~, DIS and hadron collisions; 

• LHPD and colour preconfinement; 

• heavy quark distribution and radiation; 

• jet-shape distribution. 

All these results, together with the known leading (and next-to-leading) con- 
tributions allows one to formulate multi-parton distributions in terms of a 
coherent-branching algorithms. This is the bases of most of the present Monte 
Carlo simulations for hard QCD (and beyond QCD) processes, see Q. 

An other area in which soft physics is involved is the small- x region in 
DIS and hadron collisions 13. Here soft "exchanged" gluon are involved rather 
than the soft emitted ones, as in the previous item. Enormous progresses 
have been made recently. These results are however less constrained by PT 
analysis since NP effects enter at all hard scales □. Further progresses in this 
area definitively requires further developments in our QCD understanding 
beyond PT expansions. 



3. Exact matrix element results. -Large technical progresses haizp 
been made both in analytical calculations □ and in numerical programs tj. 
These results are essential in order to make quantitative predictions away 
from the soft or collinear regions of phase space. This is achieved by matching 
these exact fixed order i^ults with the ones from resummations described in 
the previous two items til. Similar matching procedure are now developed on 
a systematically bases in Monte Carlo simulations L3. 

In the following I will illustrate the QCD calculations of jet-shape observ- 
ables which are collinear and infrared safe (CIS). They are expressed as linear 
combinations of the momenta of emitted hadrons so that their values do not 
change if two momenta are replaced by their sum when the two are collinear or 
when one is soft. For these CIS distributions one reaches the highest accuracy 
available in PT calculations. This is due to the fact that the cutoff of collinear 
and infrared singularities cancels and all PT terms are finite. Near the phase 
space boundary the cancellation are partial and logarithmically enhanced con- 
tributions enter the PT terms which needs to be resummed. The particular 
interest of these distributions is not only that they can be computed at the 
highest PT accuracy, but that they have large 1/Q-power corrections which 
are of NP origin. These corrections are phenomenologically important. Once 
we have a reliable PT estimate, with the help of data we can then constrain 
the size and hopefully the nature of these NP corrections. 

After a brief review of the results for event shapes in e + e~ annihilation, I 
discuss how these results can be extended to jet shape distributions in hadron- 
hadron and lepton-hadron collisions. This will open the way to the study of 
further multi-jet observables. 

2 Hadron emission in e + e annihilation 

Two-jet event shapes (such as thrust T, broadening B, heavy-jet mass Mh 
and C-parameter) have been intensively studied both theoretically and ex- 
perimentally. The state-of-the-art level of PT analysis of two-jet observables 
consists of Ea: resummation of all double- (DL) and single-logarithmic (SL) en- 
hanced contributions (due to multiple radiation of soft and collinear partons); 
matching of the resummed result with the exact fixed order result. To_reach 
such an accuracy one needs: soft matrix elements at two-loop order □ and 
reconstruction of the running coupling at the proper scale; factorization and 
exponentiation of observable and kinematical constraints; proper treatment 
of hard intra-jet and soft inter-jet radiation. ._. 

Besides the pure PT contributions, experimental data (see for instance 113) 
have revealed the presence of power suppressed 1/Q corrections {Q is the 



e + e~ center-of-mass energy), which emerge as a siuft both to mean values 
and differential distributions. It is widely believed til that such contributions 
are of NP origin and arise from the running of the coupling into the infrared 
domain. 

A systematic way to deal with power corrections is provided by the disper- 
sive approach lij, in which a prescription is given to extend the QCD coupling 
into the confinement region. In terms of this dispersive coupling, the typi- 
cal large distance contribution to an event shape V can be factorized in an 
observable dependent coefficient cy and a universal NP parameter ao- More 
precisely, in any 'linear' shape variable V the dependence on the rapidity rji 
and transverse momentum kt% of emitted hadron i factorizes, 



SV = Cf cy A , cy = f dr) f v (r)) , 
A NP = f & (ao(Mi) + 0(a s (Q))) , a (Mz) = £ C dk a s (k) . (2) 



Here M is a known coefficient, the Milan factor which takes into account 
the non complete inclusiveness of the observable. The NP parameter a is 
the average of the running coupling in the region k < /ij and measures the 
interaction strength in the confinement region. The 0(a s (Q)) piece in ^ is 
needed in order to merge PT and NP results in a renormalon free manner 
and ensures that the final answer is independent of the infra-red matching 
scale hi. The NP parameter ao(2GeV) has been measured and appears to-itc 
consistent with the universality hypothesis within a 10-20% accuracy, see t3. 

Only very recently these techniques have been extended to the case of 
three-jet observables in the near-to-planar region.. In particular the mean 
values and distributions of the thrust minor T m £3 and the D-parameter 113 
have been computed. Both of them arc a measure of QCD radiation out of 
the event plane. This study has revealed an unexpectedly rich geometry and 
colour dependence both of the PT result and of the NP corrections. The 
crucial point is the rapidity dependence of the observables. Actually, T m is 
independent of rapidity, while D is exponentially damped. On the PT side, 
this implies that hard parton recoil affects the observable in the T m case, 
while it does not contribute to D. Therefore the T m distribution is sensitive 
not only to the underlying hard event geometry (the angles between the jets), 
but also on its colour configuration through the kinematical constraints which 
define the event plane. 




(1) 



so that the NP contribution to V becomes 



As far as NP corrections are concerned, one finds that the shift to the 
D distribution is simply geometry dependent, while the T m distribution is 
also squeezed, since the shift depends logarithmically on T m . This behaviour 
(already encountered for the B in two-jet events) results from a complicated 
interplay between PT and NP effects. 

The same physical features appearing in the NP shift to the T m distri- 
bution are also present in hadron-hadron and lepton-hadron collisions, as we 
are going to discuss in the following. 

3 Out-of-plane radiation in hadron collisions 

We consider hadron collisions and study an observable similar to T m E3. We 
select events in which a Z$ is produced with large transverse momentum. We 
can then define an 'event plane' as the one containing the Zq momentum and 
the beam axis and introduce a measure of out-of-plane radiation 

^out = $>ri. (3) 

h 

Here p°^ is the out-of-plane momentum of the hadron h and the sum extends 
to all hadrons with rapidity in the range \rjh\ < t/q, in order to avoid measure- 
ments in the beam region. At Born level one has two incoming partons p\ 
and P2 and an additional hard parton p^ recoiling against the vector boson q. 
The event-plane definition gives rise to the conservation law 

P3 Ut + E fc ° Ut = 0, (4) 

i 

so that only p% can take an out-of-plane recoil, while the remaining hard 
momenta are fixed in the event plane. 

The main difference between hadronic and e + e~ collisions is the presence 
of initial state radiation. Due to coherence of QCD radiation, its contribution 
can be factorized giving rise to the standard parton density function for the 
incoming partons P2,P3- However, while in the total cross section the hard 
scale is the Zq hardness, here the scale is K out . The PT distribution is essen- 
tially given by the product of the initial state parton distributions -Pm (-Kout) 
and a CIS "radiation factor" 

S(X ou t) ~ P in (K out ) ■ e- R ^t) . s(R') . (5) 

The radiator R is the same as that of T m in e + e~, it resums all DL contri- 
butions and accounts for SL effects due to hard intra-jet and soft inter-jet 



radiation. It is given in terms of three hard parton 'antennae', each one pro- 
portional to the colour charge C a of emitting parton p a (C a equals Cp for a 
quark and Ca for a gluon) 

n(is \ ^ n f Qa dk2a s (2k) Q a ^ a s 2 Q a 

a=l J K out a out 

The hard scales Q a are determined essentially by soft radiation at large angles, 
which is a characteristic of multi-jet observables. In particular, apart from a 
factor due to hard collinear splitting, the scale for the quarks is the invariant 
mass of the qq system, for the gluon is its invariant transverse momentum 
with respect to the qq pair. The SL function S depends on the logarithmic 
derivative of the radiator R' and accounts for all effects coming from multiple 
secondary radiation. 

The NP correction is proportional to the same parameter which 
enters the power corrections to e + e~ event shapes in || and is given by 

SK out = ^A NP f Cifa) - m) + c 2 (vo + m) + c 3 In , (7) 

where Qt is the transverse momentum of the emitted boson. Such a shift 
arises from integration over the rapidity of emitted gluons (see eq. |^). Since 
-Kout is independent of rapidity, one has to carefully consider the effective 
rapidity cutoff. For radiation from the incoming partons p± and p 2 it is given 
by the distance between 773 (the rapidity of p^) and the experimental resolution 
r/Q. On the contrary, for a NP gluon emitted from p 3 , it is the hard parton 
recoil momentum p^ nt which provides the needed cutoff. Moreover, since p^ 
always takes recoil, one has p% ut ~ K out . This interplay between PT and NP 
emissions makes the shift logarithmically dependent on K OU f 

Unfortunately, this is not the end of the story. In hadronic collisions 
one has to add a soft contribution due to the beam remnant interactions. 
This is the so-called 'soft underlying event', which was systematically studied 
and introduced for the first time in the analysis of the 'pedestal height' in 
hadronic jet production Ed. Therefore the analysis of K out is also important 
to understand the physics of soft collisions. 

4 Out-of-plane radiation DIS 

Even in DIS one can define an event plane and measure the out-of-plane 
radiation. In the Breit frame, we define the thrust major Tm in analogy with 



e + e annihilation 

T M Q = max }\ph ■ n M \ , (8) 

n M -n=Q ' 
h 

where n and Q 2 = —q 2 are the Breit axis and the virtuality of the exchanged 
boson q respectively. As in the previous case, the sum extends to all hadrons 
not in the beam direction (with 77 < 770 ) ■ At Born level, we have one incoming 
parton p\ which is struck by a vector boson q and produces two hard large 
angle partons p 2 and p 3 recoiling one against each other. We define K out 
again as the cumulative out-of-event-plane momentum, where the event plane 
is the one containing ft and the Tj\/ axis. The event-plane definition implies 
the following kincmatical constraint 

?r +E fc ° ut =P3 ut +E fc ° u \ (9) 

with U (D) the region containing parton p 2 (^3). 

The PT -Kout distribution has the same form as ||, with Pj„ c [K out ) the 
structure function of the incoming parton, again at the scale -Kout- The radi- 
ator turns out to be the same as in e + e~ and in hadron-hadron collisions, due 
to universality of soft and collinear radiation. However the exact functional 
form of the SL function S differs between these, due to different event-plane 
kinematics. 

The NP shift is similar to the one in [7] 

9 „ n / nNP p ,„ n NP r>NP\ 

«*- - ! aNP ( c > " V- + ft ■» w\ + a in m ■ m 

with Q^ P the proper NP hard scales, which are proportional to the ones 
entering the PT radiator. The crucial difference is that here both p 2 and p% 
are not fixed in the event plane. This implies that the way the NP correction 
in [h] affects the PT distribution is very different according to which phase 
space region is considered. Namely one has the two regimes 

• a s In 2 K out /Q ^ 1: one has well developed QCD radiation, so that all 
p a (a = 2, 3) take recoil and p° ut ~ K out . This gives rise to a logarithmi- 



cally enhanced shift 



nNP e r,„ NP NP 



SK out ~C 1 \n^-— + C 2 ln^—+C 3 \n^—. (11) 

H ^out -ft out 

a s In 2 KTout/Q <C 1: radiation from one hard parton dominates. Accord- 
ing to the event-plane kinematics in|p, when PT radiation comes from the 



U region (which happens with probability [\C\ + C^)/ (2Cf + Ca)), only 
P2 takes recoil, so that that p™* ~ ^outj while p^ nt <C P2 Ut - The part of 
the shift proportional to C2 gets logarithmically enhanced as in the pre- 
vious case, while the one proportional to C3 gives rise to a very singular 
contribution proportional to 1/ \ful, coming from the average of In Q /p^ 
over its DL Sudakov form factor exp {-(|Ci + C 3 )^ In 2 Q/p^ 1 }- This 
gives 

G NP Vo nNP n 

5K out ~ d In ^— + C 2 In :g— + C 3 . (12) 

Q ^out 2^(±C\+C 3 )a s 

The converse argument holds for an emission in the D region. 

In conclusion, we have now a promising theoretical method to deal with multi- 
jet observables. Not only are we able to identify all sources of SL contributions, 
but we can also relate the NP 1/Q power-suppressed corrections with the ones 
which affect e + e~ two-jet shapes. A next step will be the extension of the 
above results to event shapes involving any number of jets. We hope that 
this analysis will greatly improve the understanding of confinement physics, 
especially at hadron colliders. 
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